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Aim 



To understand the effect of Bose-Einstein statistics on the 
thermodynamic properties of bosons. 

Objectives 

1. To derive the heat capacity of a solid using the phonon 
model. 

2. To derive the frequency spectrum of the black body 
radiation. 

3. To derive the condensation temperature of a Bose-Einstein 
condensate. 
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6.1 Phonons 
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Specific heat capacities of solids 



Previously, we have looked at Einstein's model of the 3-D 
simple harmonic oscillator, in which each atom in a solid is in 
its own potential well. The predicted heat capacity was 
accurate at high temperature, but fails at low temperature. 




In 1912, Peter Debye, a Dutch physicist working in Germany, 
produced a theory which turns out to be very accurate both 
hot and cold. 
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Atoms vibrating in a solid 



In Debye's model, we think of the atoms in a solid as 
connected by springs. They vibrate in 3D in a complicated way. 

The problem can be simplified by treating the vibration as a 
superposition of waves of different frequencies. 




If we suppose that the atoms are fixed at the edge of the solid, 
the frequency would be quantised. A quantum of this vibration 
energy is called a phonon. 
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Quantised thermal waves 



We can count the number of states for the phonons, like we 
did for atoms in an ideal gas, or electrons in a metal. 

Unlike electrons, however, the exclusion principle does not 
apply to phonons. 

So each energy state can be occupied by any number of 
phonons. 

This means that we need yet another way to find the 
macrostate. This was worked out by Bose and Einstein in the 
1920s. 
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Bose-Einsten Statistics 



The method is similar to what we have used for the ideal gas. 



There, we have made the assumption that there are many 
more energy levels than there are particles. This has meant 
that it is unlikely for two gas atoms to occupy the same state. 



We do not make this assumption now, since we would be 
particularly interested in the low temperature behaviour. 
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Since the exclusion principle does not apply to phonons, more 
than one phonons can occupy the same state. 
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Bose-Einsten Statistics 



Suppose, in the energy bundle i, that there are gi levels and 
phonons. The phonons are now completely free to arrange 
themselves among the energy levels. 

This means that we have ni + gi objects altogether in the 
bundle. There are ni of one type, and gi of the other type. The 
number of possible arrangements would be 



The total number of arrangements for all bundles would be 
obtained by multiplication of the number for every bundle: 



Tij\gj\ 



Q = n 



fii ! gi ! 
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Bose-Einstein Distribution 



Then we need to maximise InQ using the Lagrange multiplier 
method. 

Using the same constraints as before on the number of 
particles TV and the total energy U, the method would give the 
following answer: 

ni 1 

Qi exp(-Ai - X 2 €i) - 1 

where Ai and A2 are the Lagrange multipliers. 

This is called the Bose-Einstein distribution function. It looks 
almost the same as the Fermi-Dirac distribution function, 
except that it has -1 in the deniminator instead of +1. 

When we apply this to phonon, however, we need to make 
some changes. 
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Bose-Einstein Distribution 



Phonons are not real particles like atoms or electrons. The 
number of phonons is related to the vibration energy. The 
number would increase when the temperature increases. At 0 
Kelvin, there is no vibration, so there would be no phonons. 



Because of this, there would be no constraint on the number N. 
We really only have one constraint - that of constant energy U. 

The Lagrange function is InQ + At/. TV does not appear. 
Maximising this function, we get 

Hi 1 



Qi exp(-A^) - 1 
With the help of entropy as before, we would find 

1 

A = . 

k B T 
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Phonon occupation number 



Maximising Lagrange function, we have obtained 



Qi exp(si/k B T) - 1 

Following the same notations that we have used for ideal gas 
and electrons, we rewrite: 

rii as n(s)ds and 
gi as g(e)de. 

This gives 

g{e)de 



n(e)de = 



exp(s/k B T) - 1 
Next, we need a formula for the density of states g(s). 

Note that although we have previously used the same symbol 
g(s) for the particle in a 3-D box, the phonon density of states 
would have a different formula. 
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density of states for phonons 



The vibrations in a solid is described by the wave equation 

dx 2 dy 2 dz 2 v 2 dt 2 
where £ is related to the displacement of the atoms. 

http: //en.wikipedia. org/wiki/Wave_equation 



The solution has a similar form to the wave function for the 
particle in a box: 

where uo is the angular frequency. 



Substituting the solution for £ into the wave equation above, 
we would get 

i j 2 



k T ~\~ k lt ~\~ k z — 



J x ' ,v y i ,v z 9 
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density of states for phonons 



Substituting into the wave equation would be easier if we work 
out each of the partial derivatives first. For example, 
differentiating twice with respect to x, we would get 

— — = e lLJl (— k^s\n k x x) sin k y y sin k z z = —k££. 

Likewise with respect to y and z. Differentiating twice with 
respect to t, we would get 

^-f = (-uj 2 e lLjt ) sin k x xs\r\ k y ys\r\ k z z = —u) 2 ^. 
dt z 

Now we can see that substituting into the wave equation gives 

1 

v- 

which is 

.2 



k T ~\~ k lt ~\~ k~ — 



x I y • ,v z O 
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density of states for phonons 



Defining k as the magnitude of the wave vector, we have 



.2 — 7^2 i u2 i u2 

and 



k — ^ 



k 2 - — 

v z 

The boundary condition requires £ to be zero at the walls. 
Therefore, in the same way as the particle in the box, the 
wavevector would be quantised: 

u — 71x71 u — n y 7T h u — Uz7T 

rZx — i i^y — a nO Kz — • 

a J a a 

Substituting these into 

.2 7^2 _|_ 7^2 _|_ z~2 

gives 



k — ^x ^y ^ ' 



k 2 = ^(n^ + n 2 + n 2 ) 
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density of states for phonons 



Rearranging, 

k 2 a 2 



2 2 2 
— Ti T ~\~ Tin. ~\~ Ti z . 



^ 2 ,v x " 

We have also obtained above the dispersion relation 

k = -. 

V 

Substituting, 

2 2 

V a _ 2 . 2 . 2 

We may now follow the same method that we have used for 
ideal gas to count the states. Comparing with the equation for 
a sphere 

R 2 = n 2 + riy + n 2 

the total number of states with frequency less than uj is given 
by the octant of the sphere: 

7T (loo\^ 



8 V 3 J 6 \vttJ 
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density of states for phonons 



Differentiating the number of states, 

G(co) = - — 
6 \v7tJ 

we get the density of states: 
where V is the volume a 3 . 
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About polarisation 



Actually, we are not quite there yet. 



There reason is that each state correspond to a wave, and a 
wave in a solid can have 2 transverse polarisations, and 1 
longitudinal polarisation. 



Previously, we have only counted the states in terms of the 
possible wavevector (n x ,n y ,n z ). In the case of solids, the 3 
different directions of vibrations are in fact different states with 
the same wavevector (n x ,n y ,n z ). 

So we have to multiply by 3 to get the correct density of state: 

V(jj 2 duj 

g(uj)duj = 3 X 
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The Debye Frequency 



There is one more practical point to note. In a solid, there is a 
limit to the highest phonon frequency. 

This is related to the fact that the wavelength cannot possibly 
be shorter than the distance between atoms. 

Let this highest frequency be uj d . 

In a solid with TV identical atoms, there are 37V energy states, or 
normal modes. This is when we count every possible frequency, 
and every possible "direction of vibration" for each frequency. 

This can be proven mathematically, but we shall skip that. We 
can find uj d by integrating the density of states: 
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The Debye Frequency 



This equation 

>w d 3Vuj 2 duj 



Jo t>M d » = J 0 



= 3N 



2 7r 2 i; 3 

can then be solved to find the Debye frequency. Integrating 
gives 

D =3N 



2 7r 2 i; 3 

Solving for uj d gives the Debye frequency: 



f 6N7r 2 v 3 ^ 1/3 
V V 
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The phonon distribution 



As we have shown earlier, the distribution of phonons is given 
by 

f x , g(e)de 

n{e)de — ; — r 

exp(s/k B T) - 1 

We shall express this in terms of frequency uo as well. The 
phonon energy e is given by 

e = huj. 

This is a result of the phonon mode 

behaving as a simple harmonic oscillator. The phonon mode is 
simply vibration of the atoms co-ordinated in a particular way. 
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The phonon distribution 



So the energy must take the form 

E n = (n + Huj. 

Each step or quantum is hv, which would be interpreted as a 
phonon particle in our statistical treatment. The relation 

£ = huu 

would be justified by the fact that its predictions agree with 
experiments. 

In terms of the frequency u>, the phonon distribution is then 
given by 

g((jj)doj 



n(uj)d(jj = 



exp(hw/k B T) - 1 
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Debye's theory for Heat capacity of solid 



The phonon distribution is given by 



n((jj)duj = 



exp(hw/k B T) - 1 

We can multiply by ftw to get the energy at this frequency 
interval: 

fyjjn{uj)auj = . 

exp(hu;/k B T) - 1 

The total energy is then 



U — I hujn(uj)duj = 

Jo Jo 



exv(hw/k B T) - 1 
We have found earlier that the density of states is: 

g ( u ) = 3 x 

Substituting into the integral for the total energy above 
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Debye's theory for Heat capacity of solid 



We get: 

^ _ 3Vh r^D uj 3 duj 
~ 2tt 2 v 3 Jo exp(huj/k B T) - 1 

At low temperature, it can be shown that this is proportional to 
T 4 . This would imply that the heat capacity C is proportional 
to T 3 . (This has been mentioned in the topic on the specific 
heat of electrons in metal.) 

At high temperatures, the heat capacity tends to 3Nk B . 
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Is verified at high temperature 



In 1912, Debye measured the high temperature behaviour of 
copper. 




0 100 200 300 

(Beg K) 



Fig. 10' 2' 4 The molar specific heat of copper: Debye theory with Qd — 309°K 
and experimental points. (After P. Debye, Ann, Physik, vol. 2P, p. 789 (1912).) 

±$$*t in*, 



He showed that the heat capacity did tend to 3Nk B , as he had 
predicted. 
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And verified at low temperature 



In 1953, at Purdue University, Keesom and Pearlman measured 
the low temperature behaviour of potassium chloride. 



3.0* io-s 




Fig. 1. C,/r v$ for KC1. •— Our results. A— Keesom 

and Clark. 



P. H. Keesom and N. Pearlman, Physical Review, vol. 91 (1953), pp. 1354-1355 



As predicted by Debye, the heat capacity was indeed 
proportional to T 3 . 
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A note on the recommended text 



If you are using the recommended text: 

Statistical Mechanics - A Survival Guide, 
by A. M. Glazer and J. S. Wark 

please note that the Lagrange multiplier method given on page 
72 is incorrect. 

Please use the version given in these lectures. 
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6.2 Photons 
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Radiation from a very hot object 



When an object gets very hot, it can give out light. 

The temperature of a volcano lava flow can be estimated by 
observing its color. The result agrees well with the measured 
temperatures of lava flows at about 1,000 to 1,200 deg C. 




http : //en . wikipedia . org/ wiki/Black_body 
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An ideal black body 



The amount of radiation emitted depends on the nature of the 
object - its colour, whether it is smooth or rough, etc. 

An ideal black body is one that absorbs all the radiation that 
falls on it. 

No real material can do this. Soot is about the best, absorbing 
all but 3%. 

In 1859 Kirchhoff had a good idea: "a small hole in the side of 
a large box is an excellent absorber, since any radiation that 
goes through the hole bounces around inside, a lot getting 
absorbed on each bounce, and has little chance of ever getting 
out again." 

http : / / galileo . phys . Virginia . edu/ classes/ 252/black_body_radiat ion . html 

The energy of the absorbed radiation reaches equilibrium 
among different frequencies in the cavity. This results in a 
characteristic spectrum that would be emitted again through 
the hole. 
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Black Body Spectrum 




http: //en.wikipedia. org/wiki/Black_body 



As the temperature decreases, the peak of the black-body 
radiation curve moves to lower intensities and longer 
wavelengths. The black-body radiation graph is also compared 
with the classical model of Rayleigh and Jeans. 
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Photons in a box 



We can derive a formula for the black body radiation by 
considering the energies of photons in a box (like a idea black 
body). 

This follows steps very similar to the Debye theory for phonons. 

The wave equation for photons is derived from Maxwell's 
equation: 

d 2 E d 2 E d 2 E _ 1 d 2 E 

dx 2 dy 2 dz 2 c 2 dt 2 
It has the same form as the wave equation for vibration of 
atoms, except that the speed is now that of light instead of 
sound. 

So we would expect that the same formula for the density of 
states can be used. 
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Photons in a box 



Like phonons, the photon number is not fixed. We get few 
photons when the box is cold, and more photons when it is hot. 

Unlike phonons, there is no upper limit to the frequency (no 
Debye frequency). The wavelength in a solid cannot be shorter 
than the distance between atoms. The electromagnetic wave in 
a box has no such limit. 

Finally, phonons in a solid can have 3 polarisations: 2 
transverse and 1 longitudinal (like sound). Photons can only 
have 2: both transverse (e.g. light). 
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Counting photon states 



Recall that for phonons, the density of states is give by: 

Vuj 2 duj 

g(u;)du; = 3 X 

where the factor of 3 comes from the 3 polarisations of a 
phonon, 2 transverse and 1 longitudinal. 

Since photons only have 2 polarisations, both transverse, the 3 
should be replaced by a 2: 

g{u)duj = 2 x 27f2c 3 
where the sound speed ^ is replaced by light speed c. 
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Distribution of photon number 



Like phonons, each energy state can be occupied by any 
number of photons. 

So it obeys Bose-Einstein statistics. 

The number density is therefore given by the same formula 

n(e)ae = — 

exp(s/k B T) - 1 

So the number of photons in a given frequency interval is 

Vup-duj duo 
n((jj)duj = 2 x 0 Q x 



2tt 2 c 3 exp(frw/k B T) - 1 
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Distribution of photon energy 



The energy of a photon is hv. The energy in the frequency 
interval dw is then 

Voj^duo duo 
hum((jj)d(jj = 2 x hw X x 



2tt 2 c 3 exp(frw/k B T) - 1 
The energy density is given by: 

( \ * ( \ Vhoj3 1 
u(uj) = aujn(uj) = — ^-^ . 

v J y J tt 2 c 3 exp(hoj/k B T) - 1 

This is essentially Planck's law for black body radiation. 
Integrating gives the total energy: 

7T 2 Vk% . 

U = —T 

15h 3 c 3 

It can be shown that this leads to Stefan's law of radiation 

rj = aT 4 

for radiation emitted by an object at a temperature T. 
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Distribution of photon energy 



Berlin, 1899. In the Imperial Institute of Physics and 
Technology, the scientists were just trying to help the German 
lighting and heating industries. 




Otto Ernst 
Lummer Pringsheim 
1860-1925 1859-1917 



Lummer and Pringsheim measured the radiation spectrum from 
a blackbody very accurately. 
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Blackbody radiation 



The results did not look at all like what the physicists thought 
it should. 




Then Max Planck came along and derived his formula for black 
body radiation that agreed very well with these measurements. 
And the rest is history. 
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6.3 Bose-Einstein Condensation 
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Bosons 



We have seen that for phonons and photons, any number of 
such particles can occupy the same energy state. 

Such particles are called bosons. They have integer spin. 

E.g. phonon has spin 0, and photon has spin 1. 

Atoms can also have integer spin number. 

For example, ( 14 0) oxygen-14 atom and the ( 4 He) helium-4 
atom are bosons. 
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Bosons, real and virtual 



Phonons and photons can appear or disappear when an object 
is heated or cooled. They are called virtual bosons. 

Atoms, however, cannot just appear or disappear. Those with 
integer spin number are called real bosons. 

Like phonons and photons, any number of such bosonic atoms 
can occupy the same energy state. 

This gives rise to a curious state of matter called Bose-Einstein 
Condensate. 
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What is Bose Einstein Condensate (BEC)? 



When it gets cold enough all boson particles can fall into the 
ground state. 
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These particles then behave as if they are a single particle with 
a single wavefunction - all of them with the same velocity! 



PHYS393/Hock 



41 



2009-2010 



1995, USA, University of Colorado at Boulder 



Cornell and Wieman observed the BEC for the first time. 

At first, they saw the expected distribution of velocities (left) 
in a gas. Then the graph becomes a single peak (right). This 
means that all atoms have achieved the same velocity - i.e. 
become a BEC. 



ABC 




Fig. 2. False- color images display trie velocity distributor: of the cloud (A) just w'lti thermal equilibrium. The condensate fraction (moa% blue and white) is 

before the appearance of the condensate, (B) just after the appearance of the elliptical, indicative that it is a highly nonthermal: distribution, The elliptical pattern 

condensate, and (C) after further evaporation has left a $amp4e of nearly pure is in fact an image of a single, macroscopicaily occupied quantum wave func- 

condensate. The circular pattern of (he iwnoondensate fraction (mostly yellow tion. The field of view of each image is 200 pj.m by 270 iim The observed 

and green) is an indication that the vetodty distribution is isotropic, con&stent horizontal width of the condensate is broadened by the experimental resolution. 
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When do the bosons condense? 



Using the Bose-Einstein distribution that we have learnt, we can 
calculate the temperature at which bosons condense into BEC. 

Recall the expression we first derived for the Bose-Einstein 
distribution: 

n% 1 

g { ~ exp(-Ai - X 2 Si) - 1 

When applying this to phonons, Ai was left out because the 
number of phonons is not fixed. 



Now that we are dealing with atoms, which are real bosons, the 
number is fixed and we have to keep Ai. This is usually written 
in terms of \± as follows: 

g(e)de 



n(s)ds — 



exp((e-fi)/k B T)-l 
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The density of states 



We have seen a number of different density of states. 

For electrons, we need to multiply by 2 because of the spin 
states. For phonons, we multiply by 3 for the three possible 
polarisationz. For photons, we multiply by 2 for the two 
polarisation states. 

So what do we do for atoms that are bosons? 

The answer is: we use the same density of state for the ideal 
gas: 

^)=^ r (2ms) 1 / 2 

This is the very first one that we have seen, before we have to 
include the additional effects of spin and polarisation. 
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The chemical potential 



/I is called the chemical potential. It is determined by the total 
number of particles. We get this by integrating: 



Jo exp((e - v)/k B T) - 1 

I have assumed that the ground state energy is zero, which is 
why I integrated from zero. 

This is not easy to solve. Instead, we shall make use of some 
approximations at low temperature. 

For a start, note that \± could depend on temperature, since 
the above equation contains T. 
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The occupation number 



We focus on the occupation number in the integral: 



exp((e - n)/k B T) - 1 
First, note that \± cannot be positive, or else for small energy 
the exponential function would be less than one. Then the 
denominator would be negative. This means negative 
occupation number, which is not physical. 

Next, we know that when temperature is low enough, a large 
number of atoms would go into the ground state e = 0. The 
number would be as large as N, which is of the order of 1 mole 
(« 10 24 ). 

This means that /(0) would be very large, so that \i must be 
very close to zero. Then the exponential function would be 
close to 1 and the denominator would be close to zero. 
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Number of atoms in the excited state 



So at very low temperature, it is safe to assume that jjl = 0. 
We can then write 



I have added a subscript ex to N. The reason will become clear 
in a moment. First, integrate this with the help of the Table of 
Integrals: 



Notice that T is in the denominator. This suggests that the 
total number of atoms decreases with temperature - that they 
are disappearing ! 

In fact, the above integral includes only atoms above the 
ground state e = 0. The "missing" atoms are in fact going into 
the ground state. 
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The condensation temperature 



The formula for N ex is only able to compute the number of 
particles above the ground state. This is because the density of 
states g(s) in the integral is zero when e = 0. 

When a substantial fraction of the particles start going into the 
ground state, this has to be added separately. So N ex is the 
number of atoms in the excited states. Hence the subscript ex. 

We know that as soon as N ex becomes less than the original 
number TV, then (N — N ex ) atoms start going into the ground 
state. 

Therefore, condensation takes place when N ex = N. 
Substituting this and solving for T, we get 
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Number of atoms in the ground state 



Above T BE , there are very little atoms in the ground state. 

Below T BEl this number increases until all atoms fall into the 
ground state. 



N 




Glazer and Wark, Statistical mechanics, p. 101 
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How cold is a BEC? 



When Cornell and Wiemann observed the first BEC in 
Colorado, they used rubidium gas with a condensation 
temperature of 1.7 x 10 -7 K. 

This is far below the temperature that can be reached using 
liquid helium (melting point 4.2 K). 

They have to use a special method called laser cooling. 

For more information on laser cooling, see: 

http : 

//physweb . bgu . ac . il/C0URSES/LAB_C/Laser°/ 0 20Cooling 0 / 0 20and 0 / 0 20trapping/colorado 0 /,20mot . pdf 
http: //en.wikipedia. org/wiki/Laser_cooling 
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What we have learnt so far. 



1. The Bose-Einstein distribution function is given by 



exp[( £ - n)/k B T] - 1 
2. For phonons and photons, the following form is used 

/GO = 1 



exp(s/k B T) - 1 
3. The density of states can take different forms: 



phonons 



photons 



bosonic atoms 



Vuj duo 
g(uj)d(jj = 3 X 

27T Z V^ 



, _ Vuj 2 duj 

g(uj)duj = 2 x 
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What we have learnt so far. 



electrons 



classical ideal gas 



<?(» = 2 x — ^(2m£) 1 / 2 



5 (s) = ^ r (2ms) 1 / 2 



3. At low temperature, the specific heat capacity for the Debye 
model tends to T 3 . At high temperature, it tends to 3Nk B . 

4. The energy distribution in a photon gas is given by: 

( \ * r ^ VhujS 1 



tt 2 c 3 exp(hu;/k B T) - 1 

This also describes the black body radiation spectrum. 

7. The condensation temperature for a Bose-Einstein 
condensate is: 

_ h 2 / N \ 2 / 3 
TbE ~ 27rmk B V2.612W 
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6.4 Exercises 
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Exercises 



Exercise 1 

Sketch the density of states as a function of frequency for 1, 2 
and 3-D solids in the Debye approximation. Consider a 2-D 
solid. Show according to the Debye theory that the heat 
capacity varies as T 2 at low temperature. In graphite it is found 
that C v oc T 2A at low temperature. Explain this behaviour. 
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Exercises 



Recall the steps leading to the density of states of a phonon in 
a 3-D solid. First, we obtain the quantised energy states: 

2 2 

u a _ 2 . 2 . 2 

Then we compare with the equation for a sphere 

R 2 = n^ + riy + n 2 

For positive integers (n x ,n y ,n z ), the total number of states 
with frequency less than uj is given by the octant of the sphere: 

= - f-i? 3 ) = - f-f 
8 V 3 J 6 \vttJ 

Differentiating, we get the density of states: 

g(u) = 3 x 

where V is the volume a 3 . There is an extra factor of 3 
because in 3-D, there can be 3 modes for each phonon 
frequency: 2 transverse and 1 longitudinal. 
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Now we repeat the steps for 2-D. First, we obtain the 
quantised energy states (no n z )\ 



2 2 

uj a 2 f 2 



TXlrp | TXln 



Then we compare with the equation for a circle (not sphere) 

R 2 = n% + riy 

For positive integers (n x ,n y ), the total number of states with 
frequency less than uj is given by the quadrant of the circle (not 
octant): 

G(u) = -(nR 2 ) = - — 
4 4 Vi;7r/ 

Differentiating, we get the density of states: 

, . Auj 
g(uj) = 2x - — ^ 

where A is the area a 2 (not volume). There is an extra factor 
of 2 because in 2-D, there can be 2 modes for each phonon 
frequency: 1 transverse and 1 longitudinal. 
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Finally we repeat the steps for 1-D. First, we obtain the 
quantised energy states (no n y )\ 



2 2 
uj a 2 



Then we compare with the equation for a line (not circle) 

R 2 = n 2 

For positive integers n x , the total number of states with 
frequency less than uj is given by half of the line (not quadrant): 

G( W ) = \{2R) = ^ 
Differentiating, we get the density of states: 

g(uj) = 1 x — 

where L is the length a (not area). There is no additional 
factor because in 1-D, there can be only be 1 mode for each 
phonon frequency: the longitudinal. 
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Summarising, we have obtained the following results for the 
density of states: 



1- D: g(u) = constant 

2- D: ^f(cj) oc uj 

3- D: ^f(cj) oc up- 



The internal energy is given by: 

^ / W D hujg(uj)duj 



j o exp(huj/k B T) - 1 
In 2-D, 

g(oj) oc a;. 

So 

(jp-dijj 

U oc 



f o exp(hoj/k B T) - 1 
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We have obtained 

U oc 



'o exp(huj/k B T) - 1 
Let 



k B T 

At low temperature, T — > 0 and x — > oo. The integral for [/ can 
be written as 

U ocT . 

JO e x - 1 

The factor of T 3 appears because of the factor of T in the 
substitution for u\ 

k B T 

UJ — X. 

h 
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Since the definite integral here 

i3 r°° x 2 dx 



r 0 e x — 1 

is a constant, we have 

U oc T 3 . 

Differentiating, the heat capacity is 

C v oc T 2 . 

Graphite has a layered structure. Each layer consists of carbon 
atoms strongly bonded, with weak forces between layers. 
Vibration (like sound) would propagate much faster along the 
layer than across the layers. So it can be thought of as 
intermediate between 2-D and 3-D. 

This explains why the heat capacity of graphite is 

C v oc T 2A . 
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Exercise 2 

According to Debye's theory, the total energy of phonons in a 
solid is given by 

^ _ 3Vh r^D uj 3 duj 

~ 2tt 2 v 3 Jo exp(huj/k B T) - 1 ' 

Show that at low temperatures, i.e. T <c hu; D /k B , Debye's 
theory predicts that the heat capacity of a solid is proportional 
to T 3 . 



[ You are given that 



•oo x 3 7T 4 

ax = 



f 0 e x - 1 15 
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To use the formula given, let the argument of the exponential 
function be x: 



k B T 

Make frequency the subject: 

k B T 

UJ = X. 

h 

Then substitute into the integral for the energy: 

4 _3 



_ 3Vfi fk B Ty r%D x^dx 
~ 2tt 2 v 3 \ H J JO e x - 1 



where the new limit is related to the old one by the substitution 
formula: 
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The condition given is T <^hw D /k B . Rearranging, we have 



The right hand side is just x D , so 

1 <C x D . 

This means x D is large. We may approximate the integral by 
replacing the limit by infinity: 




Then using the formula given 



U = 



2tt 2 v 3 V h 



3VH fk B T 
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We have found the total energy. 

Vn 2 ktT 4 



U = 



B 



ioh 3 v 3 

Note that this is proportional to T 4 . To get the heat capacity, 
we differentiate with respect to temperature T: 

2VtT 2 /c 4 T 3 



C = 



B 



5h 3 v 3 

which is proportional to T 3 . 
This answers the question. 
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Some remarks: 



Going back to our approximation, how large must x D be for the 
approximation 

%X D x 3 dx ^ r°° x 3 dx 
lo e x - 1 ~ JO e x - 1 

to be valid? We can get some idea if we carry out the 
integration numerically: 



x D 


c oo doc 
J e x_-\_ 


error 


1 


0.225 


96% 


10 


6.432 


1% 


100 


6.494 


io- 18 


oo 


7T 4 /15 


0 



As long as x D is more than 10, the approximation would be 
better than 1%. 
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Recall that x D is related to the Debye frequency by 

XD = ^ 

This means that for 1% accuracy, we should have 



k B T 



> 10. 



Rearranging, we get 



where 



T < — 
10 



"d — 



k B 

6 D is called the Debye temperature, and can be used as a 
reference when we want to tell whether a temperature is high 
or low. Note that in the numerator is hup, which is the highest 
possible phonon energy. 
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Exercise 3 



The energy distribution of a gas of photons in a box given by: 



tt 2 c 3 exp(hw/k B T) - 1 
Wien's displacement law for black body radiation is given by 

\ _ b 

Amax — — , 

where Amax is the wavelength of the peak of the radiation 
intensity, and b is 2.90 x 10 -3 m.K. 



Derive Wien's displacement law from the energy distribution. 
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We must rewrite the distribution in terms of wavelength A: 

c 



Differentiating, 



J A 



d(jj = 7z-d\. 

A 2 



Substituting into the energy distribution: 



tt 2 c 2 A 5 exp(hc/\k B T) - 1 

To find the maximum, we can ignore the constants in u(u>) and 
just focus on the part that depends on uj: 

1 1 



A 5 exp(hc/\k B T) - 1* 
To find the maximum, differentiate this ... 
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and set to zero: 



5 



1 



1 (-hc/X 2 k B T) exp(hc/Xk B T) 



= 0. 



A 6 exp(hc/Xk B T) - 1 A 5 [exp(hc/Xk B T) - l] 2 
Multiplying by X 6 [exp(hc/Xk B T) - l] 2 , 

he 

-5[exp(hc/\k B T) - 1] + — — exp(hc/\k B T) = 

\k B l 

We can simplify this by defining 

he 



x = 



Xk B T 



Substituting, we get 



5(e 



X 



1) + xe x = 0. 
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If we plot each side of the equation, 

5(e x -l)=xe 
we see that there are 2 solutions. 



X 



1200- 



1 0DD- 



BOO' 



600- 



400- 



200- 



-200- 





y= 5"Cexp«|- 1) 




y = x .* exp(x) 




7 






1 I 1 1 I 1 1 1 1 1 I 1 1 1 


— i 1 1 1 1 



We can solve this numerically. The answers are 

x\ = 0 and x^ — 4.965. 

PHYS393/Hock 70 



2009-2010 



Exercises 



The nontrivial solution is 

x 2 = 4.965. 
Substituting this into the definition for x, 

h 

x — 



Xk B T 

or making A the subject, we get 

he 

A — — 



Comparing this with the displacement law, 

\ _ b 
Amax — — , 

we see that 

he 



b = 

Substituting the constants and calculating, we get the Wien's 
displacement constant 2.90 x 10~ 3 m.K. 
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Exercise 4 

Bose-Einstein condensation is achieved with a gas of rubidium 
atoms, which have an atomic weight of 85.47. The 
condensation temperature is 1.7 x 10 -7 K. What was the 
number density of the condensate? 

[ Atomic mass unit u is 1.6605 x 10 -27 kg.] 
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Recall the formula for the condensation temperature: 

h 2 ( N \ 3 / 2 



T BE — 



) 



2irmk B V2.612V 
We can rearrange this to get the number density: 

N ^ _ {27rmk B T BE \ 3 / 2 



= 2.612 



I *2 



V 

The condensation temperature is given as 

T BE = 1.7 x 10~ 7 K. 

The mass is given by 

m = 85A7u = 85.47 x (1.6605 x 10~ 27 ) kg. 

Substituting these into the above equation for the number 
density, we get 2.74 x 10 19 m -3 . 
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Exercise 5 



Show that below the condensation temperature, T BE , the heat 
capacity of a gas obeying Bose-Einstein statistics is given by 

T \3/2 

C v = 1.93Nk B 



BE 



[You are given that 

•oo x 3 / 2 



J 0 e x — 1 



-dx = 1.78, 
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Recall that at very low temperatures, the chemical potential n 
is approximately zero, and the number of particles above the 
ground state is given by: 

N = f 00 g(s)ds 

ex Jo exp(e/k B T)-l 
To obtain the energy, we insert e into the integral: 

u= f°° g(e)ede 

Jo exv(s/k B T) - 1 

Substituting the density of states 

/ N A-mnV . .1 /9 

gives 

_ 4m7rV f<>° V2ms 3 / 2 ds 



h 3 Jo exp(e/k B T) - 1 
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Define 

£ 

X = 



k B T 

and substitute into the integral for energy: 

h 3 * Jo e x - 1 

We can now use the given result: 

-dx = 1.78. 



! 0 e x - 1 
This gives: 

U = ^^V2^(k B T)^ x 1.78. 

The heat capacity can then be obtained by differentating this 
with respect to temperature. First, we express this in terms of 
the condensation temperature. 
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The condensation temperature is given by: 

_ h 2 f N \ 3 / 2 
TbE ~ 27rmk B V2.612W 
Make V the subject 

2.61 \27Tmk B T BE J 
and substitute into the above result for U: 

2 1.78 T 5 /2 
U = Nk B — 

Then differentiating with respect to T: 

3/2 

C V = 1.93Nk B ' 



T BE, 

we get the heat capacity for the Bose-Einstein condensate, 
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